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We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, 
and show that the frequency dependence of the long- wavelength spin response is well constrained by 
sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules 
are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two- 
particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, 
the sum-rules suggest that the strength of the spin response peaks at uj ~ 40-60 MeV, decays rapidly 
for u) >100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy 
may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of 
relevance to core-collapse supernova. 

PACS numbers: 97.60.Bw, 26.50.+X, 95.30.Cq, 26.60.-c 

I. INTRODUCTION 

The spin response of dense neutron matter plays an essential role in determining neutrino interaction rates in 
neutron stars and supernovae [TH5]. (For the effect of spin response on photon interaction with nucleon magnetic 
moments, see discussion in Ref. [B].) Since the energy and momentum transfer between neutrinos and matter is small 
compared to the Fermi energy and momentum, degeneracy and many body effects can strongly modify interaction 
rates. The spin response of neutron matter is an intriguing problem in that a non-zero response requires the coupling 
of spin and space through the tensor and spin-orbit components of the nuclear force. 

We study the response in the specific limit of zero temperature and zero momentum transfer, and discuss how this 
limiting case will be useful to understand the more general behavior encountered at finite temperatures in neutron 
stars and supernovas. At zero temperature the spin response can be obtained through a combination of sum-rules 
and a calculation of the high-energy part of the response. The sum-rules and the high-energy behavior resolve nuclear 
interactions with momenta of the order of the pion wavelength, and we use nuclear Hamiltonians previously found to 
be reliable in describing relevant excitations and their coupling to the ground state in the other contexts. 

Bremsstrahlung reactions such as n + n — > n-\-n+i>+v are an important source of neutrino pair production in dilute 
neutron matter. When neutrons are non-relativistic, neutrino emission occurs primarily due to fluctuations of the 
nucleon spin. Density and current fluctuations are suppressed by the velocity due to particle number and momentum 
conservation. This dominance of spin fluctuations is unique feature of nuclear systems because strong non-central 
tensor and spin-orbit forces that do not commute with the spin operator lead to enhanced spin fluctuations even in the 
long- wavelength limit (q — > 0). Its importance in neutrino production rates was first realized in pioneering work by 
Friman and Maxwell [1 . They calculated the neutrino production rate in the long- wavelength limit using the one-pion- 
exchange (OPE) potential in leading order perturbation theory (Born approximation). In subsequent work, Hanhart 
et al. calculated the neutrino production rate in neutron matter using a low-energy approximation that relates the 
rate directly to observed nucleon-nucleon phase shifts obviating the need to rely on either perturbation theory or a 
specific choice for the nucleon-nucleon potential [7] . While these calculations have provided a useful benchmark, they 
neglect many-body effects and their regime of validity is restricted to weakly correlated dilute systems. 

The inclusion of many-body effects have relied on diagrammatic perturbation theory where specific corrections to 
long-distance and long-time behavior of nucleon propagation in the medium are incorporated. The finite lifetime of 
quasi-particles, screening of the weak axial charge, as well as screening of nucleon-nucleon interactions due to particle- 
hole polarization effects at finite density have been investigated by several authors [5J IMT2] . These calculations have 
shown that these corrections are important and generically tend to decrease neutrino production rates. On the other 
hand, attempts to include in-medium softening of the pion propagator and corrections to the nucleon propagators 
and weak vertices j!3H15j have shown that the neutrino emissivity can be enhanced. However, these methods neglect 
terms in many-body perturbation theory and it is presently difficult to estimate associated errors. To overcome this 
shortcoming, we adopt a different strategy, and use Quantum Monte Carlo (QMC) to compute three lowest order 
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sum-rules which are described below in { III We supplement these sum-rules with asymptotic form of the two-particle 
response valid at high frequency to deduce the distribution of strength of the spin response function at lower energies 
of relevance to astrophysics. 



II. NEUTRINO EMISSIVITY AND THE SPIN STRUCTURE FUNCTION 



From the point of view of many-body theory, neutrino interaction rates in the medium can be factored into a 
product of two terms: (i) the correlation functions of the dense medium, and (ii) kinematical factors and coupling 
constants associated with neutrino currents. The latter are well-known and relatively simple functions of the neutrino 
energy and momentum. In contrast, the spin, density and current correlation functions are complex functions of 
temperature, density, and the energy and momentum transfer because multi-particle dynamics and correlations in the 
ground state of the strongly interacting system play a critical role. 

The dynamic spin structure factor S a {io, q) of neutron matter encodes the linear response of neutron matter to spin 
fluctuations and is defined as [2] 
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where s(t, q) = V^ 1 J^ =1 e~ %q ' ri ^ <Xj and cr, is the spin operator acting on the i th nucleon at time t. The second line 
expresses the same response as a sum over final states |/) coupled to the ground state through the time- independent 
spin operator. 

Alternatively in terms of the field operators, s(t, q) is Fourier transform of the spin density operator s(x) = 
~ip + (x) t ip{x) with r being the usual Pauli matrix and ip(x) is the non-relativistic field operator. The normal- 
ization factor 4/3n where n is the neutron number density ensures that the dynamic form factor is canonically 
normalized such that S(q — ¥ oo) = 1 for the non-interacting Fermi systems and conforms to the standard definitions 
of the sum-rules discussed below in §III| 

The rate of neutrino pair production can be expanded in powers of the nucleon velocity and the momentum of the 
neutrino pair [T]. The neutrino emissivity of neutron matter denoted by Q, and defined as the rate of energy loss due 
to neutrino pair production per unit volume and per unit time, to leading order in the neutron velocity and neutrino 
momentum is given by[9] 

Q = ^^fdu^e-^S^), (2.2) 

where Gf — 1.18 x 10~ n MeV~ 2 is the Fermi constant of the weak interaction, Ca = —1.26/2 is the neutron 
neutral-current axial coupling constant. Note due to the strange-quark contribution to the nucleon spin, Ca may be 
modified in neutral current processes by few percent in the energy range of interest to supernova [IB] . Here we have 
not included this modification for simplicity. At low temperature, when T -C Ep n , where E-p n is the neutron Fermi 
energy, the neutrino pair momentum q is small compared to the both the Fermi momentum kp n and the intrinsic 
momentum scales associated with the strong interaction, and may be neglected and S^^ui) — S q = 0). Hence in 



Eq. (2.2) only S a (u)) appears and it is both a function of density and temperature as implied by the ensemble average 



denoted on the RHS of the equation. 



III. SUM-RULES 



The spin response describes the coupling to the ensemble of final states obtained by flipping all the ground-state 
spins in neutron matter. If spin and space are uncoupled, spin is a good quantum number and there would be no 
response at zero momentum transfer. However, the spin-orbit and tensor interactions (acting only in relative p— waves 
and higher in neutron matter) induce a finite expectation value of (S 2 ) even at T=0 and a finite response results. 
The spin-orbit and tensor interactions are of pion range or less, so they predominantly affect neutrons coupled to spin 
1 at a pair separation typical for nearest neighbors at that density. Although there is zero total momentum transfer, 
the two interacting particles can nevertheless have significant relative momenta in the relevant final states. 

The overall strength and energy distribution of the response can be characterized through the relevant sum-rules. 
We employ QMC to compute the low order sum-rules that relate moments of S„(u>,q) to its ground state properties. 
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We then combine these sum-rule constraints with asymptotic high-energy behavior expected in the two-particle system 
to obtain constraints on the distribution of strength of S a (uj) as a function of uj at q = 0. For the same reason, the 



response in Eq. (2.1) is solely due to the excitation of multi-particle states as single particle excitations vanish for 
these kinematics . 

Though we ultimately desire information about the spectrum and coupling to the excited states of the system, the 
moments of the sum-rules defined by the relation 

/>oo 

S'2 = S a {u,q = 0) u n duj, (3.1) 
Jo 

are calculable as ground state properties. The sum-rules provide a simple and systematic means to eliminate explicit 
dependence on the intermediate excited states of the system. The relevant excited state information is sampled by 
operators contained in the nuclear Hamiltonian. In this study we use the following sum-rule relations: 

5- = | (3.2) 
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where \<y — dn^/d/ia is the spin susceptibility of the interacting ground state |0) of the nuclear Hamiltonian iJ/v, and 
n a and /v are number density and chemical potential of particles with spin a (±1/2). Our strategy here is to evaluate 



the right hand side of Eqs. (3.2), (3.3) and (3.4) using QMC and use this information to constrain the behavior of 
S(lu) for values of uj relevant to the calculation of neutrino production. 

This strategy is not new, in Ref. [8] estimates of the and sum-rules were used to argue that spin response 
function must saturate at high density, and in Ref. [17] . sum-rules were used to estimate the relative importance of 
multi-particle excitations to the response function in the kinematical regime where uj > q. Our work improves upon 
these earlier studies in two respects: (i) we compute and combine for the first time all three sum-rules to constrain 
both low-frequency and high-frequency behavior of S(uj,q = 0); and (ii) we deduce the high-frequency response or 
short-time behavior of the two-particle dynamics where they dominate in the many-body system by direct calculation 
of the two-particle matrix elements. 

To compute the expectation values of operators in the ground state needed to evaluate the sum-rules we use a 
non-relativistic nuclear Hamiltonian with local 2-body potentials of the form 
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where Ofj 1,4 = (1, <Ti-(Tj, Sij, L- S), and Sij — (3cri ■ f <x 2 -f— ovo^) is the tensor operator, and L- S is the spin-orbit 
operator. We employ the Auxiliary Field Diffusion Monte Carlo (AFDMC) method described in Ref. [HI [19] and use 
the Argonne AV8' form for the two-body interaction as it provides a good description of properties of light nuclei |20j . 
The AFDMC calculations use auxiliary field quantum Monte Carlo techniques to treat the spin and spatial degrees 
of freedom in neutron matter. They have been used extensively to calculate the equation of state of neutron matter, 



and also the spin susceptibility [21]. We use AFDMC to compute the sum-rules expressed in Eq. (3.2 3.3) and (3.4). 
Note the static structure function S° and energy-weighted sum rule have been previously evaluated for Argonne 
potentials [22 , first based on variational methods [551 HI] • 

The S~ l sum-rule is calculated by considering the ground state of neutron matter in the presence of a magnetic 
field as proposed in Ref. [21] , The energy of neutron matter in the presence of a magnetic field is: 



E(p) =E(0)-bP+{l/2)P 2 E"(0), (3.6) 

where E(0) is the ground state energy in the absence of a magnetic field, P = is the spin polarization, and 

the spin susceptibility x<j is 

X, - fP^y (3-7) 

The calculations are performed for zero magnetic field and a finite magnetic field for of order 60 particles in periodic 
boundary conditions. The system we simulate has finite number of up and down neutrons, and the magnetic field is 
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chosen in such a way the finite system is close to the thermodynamic limit as described in Ref . |21j . In a non-superfluid 
system, the calculation of the spin susceptibility yields the S~ x sum-rule. 

We calculate S® by computing the spin-dependent pair correlation function and evaluating the structure function 
at q = 0. The spin correlation function is defined by 

where ip is the ground state of the system. The AFDMC method is useful to compute the expectation values of mixed 
operators like (^t\0\i/j). We use Variational Monte Carlo (VMC) to extrapolate the value of operators that are given 
by (O) = 2(0) m i X — (0) vmc as described in Ref. [THUS]. The resulting g(r) is used to obtain the structure function 
S°(q). We show g a (r) and S°(q) in Fig. [j] We finally evaluate S° sum-rule by taking the q — > limit as indicated in 
Eq. @. 

The energy weighted-sum-rule can be calculated by the expectation value of the tensor and spin-orbit interactions 



when q = 0. For the Hamiltonian of Eq. (3.5 ) we have 



S* 1 = -3^ E PMrMi) + (Mrv)L ■ S)) • (3-9) 

i<3 

Because the variational wave function ^>t used as input for AFDMC contains neither tensor nor spin-orbit corre- 
lations, the most accurate way to obtain these expectation values is by calculating the energy as a function of the 
spin-orbit and tensor interaction strengths and using the slope of the energy with respect to these couplings to produce 
the true ground-state expectation values. 

These initial calculations are performed with the AV8' NN interaction without any three-nucleon interaction. Based 
upon simple estimates of the strength of the three-nucleon force, we would expect of order 10 — 20% corrections to 
the sum-rules from the three-nucleon interaction. We are exploring this dependence and will report these results 
separately. 

In computing the ground state properties in AFDMC we neglect the role of pairing and superfluidity. This will 
restrict our study to the calculation of the neutrino emissivity at temperatures that are large compared to the neutron 
pairing gaps in neutron matter but still small compared to the Fermi energy. Thus, our results will be applicable to 
ambient conditions in the supernova but will not apply to old neutron stars where neutron matter is likely to be below 
the supcrfluid critical temperature. For T <C A where A ss 1 MeV is the superfluid gap, the number of quasi-particles 
is exponentially suppressed and response is vanishingly small. In vicinity of the critical temperature, Cooper pair 
breaking and formation, as well as collective modes can enhance spin-fluctuations at a frequency u; ~ (1 — 2)A 26J. 
The spin response function and the neutrino emissivity in the superfluid phase is expected to be qualitatively different 
and is dominated by the pair recombination processes and the decay of finite energy collective modes |27l 128) . It may 
be possible in the future to examine this regime more critically using techniques similar to those developed here. 

The AFDMC results for the sum-rules are shown in Table |TJ where the individual sum-rules and average excitation 
energies defined by u>o — S^/S^ 1 and Q\ = S^/S® are listed. The density dependence of the 5° sum-rule is quite 
modest over the range of densities considered. 



Table I: AFDMC results for the sum-rules 



Density (fm 


S- 1 (MeV- 1 ) 




(MeV) 


wo (MeV) 


Oi (MeV) 


n = 0.12 


0.0057(9) 


0.20(1) 


8(1) 


35(9) 


40(8) 


n = 0.16 


0.0044(7) 


0.20(1) 


11(1) 


46(11) 


55(8) 


n = 0.20 


0.0038(6) 


0.18(1) 


14(1) 


47(12) 


78(10) 



The spin susceptibilities shown in table [I] correspond to x/xf = 0.37,0.34, and 0.34 for p = 0.12,0.16, and 0.20 
fm~ 3 , where xf = rnkp/T^ 2 is the spin susceptibility for free fermi gas. At the lowest density this is very similar to 
results obtained in [5T] , at the highest density our result is approximately 20 per cent lower for the susceptibility. The 
difference may lie in the fact that the three-nucleon force used in |21j is repulsive in unpolarized neutron matter, and 
less so in spin-polarized matter. 

The average energies Cu and u>i are extracted from the sum-rules as estimates for the energy of the peak of the 
response, and their difference is a measure of the width of the distribution. The fact that the calculated u>o and uj\ 
are fairly similar indicates a moderately narrow peak in the response. A positive definite response requires u}\ > loq. 
The peaks shift to higher energy with increasing density, as expected. The tables also indicate that the strength 
distribution gets more diffuse with increasing density with strength being pushed out to higher energy. 



5 




Figure 1: (Color online) The static structure function S®(q) computed at saturation density. In the inset we show the 
corresponding spin pair correlation function ga(r). Free particle results are given as dashed lines. 



IV. ASYMPTOTIC FORM AT HIGH ENERGY 



In order to constrain the low-energy response relevant for astrophysical applications using sum-rules we need some 
knowledge of the behavior of S a (u) at large u>. In this regime the response probes the short time behavior of the 
many-body correlation function and on general grounds we can expect this to be dominated by two-particle dynamics. 
This intuitive expectation can be cast in more formal terms using the operator product expansion originally developed 
by Wilson as a standard technique in quantum field theory. The operator product expansion has been used to analyze 
short-time behavior of the density- density correlation function in a strongly interacting non-relativistic fermi gas 
[291 130) . Adapting this to the spin-density operator, the relevant expansion in this case organizes S a {uj) in terms of 
local operators in inverse powers w, and is given by 

J dt e wt J d 3 x ^<rip(t,R+x) ^atp(0,R-x) = tWi(w) {1) {R) + iW 2 (uj) {2) (R) + --- (4.1) 

where the expectation value of the local operators 0^ n \R) depends on the many-body ground state but the Wilson 
coefficients Wi(u>) depend only on few-body physics with i incoming and outgoing asymptotic states. For q = the 
Wilson coefficient Wi(w) vanishes identically in spin saturated system and the leading contribution is due to 
The functional form of W^iiiS) is determined by the matrix elements of the spin operator between two-body scattering 
states. This implies that up to an over all constant which depends only on the ground state, S a {u>) at high frequency 
is determined by the the two-body matrix elements. In general, this will depend sensitively on the short-distance 
behavior of the two-nucleon interaction and will be model dependent. However, to extract the response at low energy 
in a model independent fashion it suffices to use in the two-body calculation, the same Hamiltonian employed in the 
calculation of the sum-rules in many-body calculation. 
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The spin response function S(q, uj) for two neutrons are evaluated as follows, 

S(q,w) = \<iP F \6 A \^i>\ 2 6(Lo + E I -E F ). (4.2) 

For spin response at g=0, the operator is the sum of spins, Oa = <?i + oi- ipi and tpp are the eigenstates of two 
neutrons in spin-triplet states and take ipi to be the ground state. 

We have calculated these matrix elements using the same nuclear Hamiltonian employed in the AFDMC by solving 
the Schrddinger equation for two-neutrons with simple box boundary condition. These results indicate that the high 
frequency behavior denoted as S£ lgh (w) is determined by two-body physics and has the following asymptotic behavior 

S^(u,)~(^)\ (4.3) 

where the density dependent quantity ui c ~ 100 — 150 MeV for the range of densities considered here and for the 
nuclear interaction used we find that i«9. 

As mentioned earlier the high frequency response will depend on model for nucleon-nucleon interactions at short- 
distance. For a correct description of the response at uj > 100 MeV, the inclusion of two-body currents and explicit 
pion and A— isobar degrees of freedom is likely to become important. However, since they are absent both in the 
many-body and two-body calculation, their consistent omission ensures that we can still obtain useful constraints on 
Sv(oj) at lower ui of interest without these ingredients. 

Using the two-body axial currents adjusted to reproduce measured tritium /3 decay |31j . we calculated the contri- 
butions to the static spin sum rule of Eq. ( |4.2| ) at q = due to the most important two-body currents - the axial 
7r-exchange A-excitation current and 7r-exchange (pair) current. It was found to be a few percent of the total static 
spin sum rule. Therefore we expect the contribution of two body currents to the dynamic spin response function at 
zero momentum transfer to be around a few percent as well. 



V. LOW ENERGY FORMS FOR THE RESPONSE 



In the regime where neutron matter behaves like a Fermi liquid, the low-energy form of the response should be 
describable in terms of quasi-particles, though the coupling of the ground state to the quasiparticle pairs as well as 
the quasiparticle interactions may renormalize quantities in the calculated response. At q — 0, a low-frequency form 
for S a (u) has been computed in Refs.[3l[5] using the quasi-particle approximation and is given by 



N(0) LOT a 

nir (1 + G ) 2 + (ujr a y 



(5.1) 



where the frequency dependent relaxation time T a {uo) is the time-scale for damping of spin fluctuations, and iV(0) 
is the density of states at the Fermi surface and Go is the Landau parameter that encodes the spin susceptibility 
mentioned earlier Eq. (3.2 1. This form of the response incorporates collisional broadening and mean field effects but 
it is mostly sensitive to the T a (uj). Note this type of functional form incorporates higher order terms in the scattering 
and thereby takes into account the Landau-Pomeranchuk-Migdal effect [3"2Tl35| . 

The spin relaxation time t ct is related to the quasi-particle scattering amplitude A C r 1 ,a- 2 (k, k') and is given by [5] 



2,3,4 



r a (uj + ei) 
where \A\ 2 = — ^ Tr 



2tt \A\ 2 T S(uj + ei +e 2 -e 3 - e 4 )5{P 1 + P 2 - P 3 - P 4 ) 



An,^ (k, k')<xi [(o-l + <r*f , Ar 1>t r 2 (-k, k')] 



(5.2) 
(5.3) 



where the first sum is over the momenta of all initial states of particle 2 and all final states of particle 3 and 4 and 
T = — /aXl — /*) + (1 — /2)/a/4 Pauli blocking factors where fa is the Fermi-Dirac distribution for particle i in the 
reaction 1 + 2 — > 3 + 4 with incoming momenta Pi and P 2 an d outgoing momenta _P 3 and P4 and relative momenta 
k = Pi — i-3, k' = Pi — P4. The squared matrix element in Eq. (5.3) is a sum over the spin projections j = 1,2,3 
of the spin operators <Ji and <T2 acting on nucleons 1 and 2 respectively. In the limit of ui <C Ep and T <C Ep, the 
appropriate average relaxation time from Eq. (5.2) reduces to 



1 



(5.4) 
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where C a characterizes the strength of non-central interactions and depends on the ambient density and T is the 
temperature. C a has been calculated using different models for the nucleon-nucleon potential for a range of densities 
in Ref. [33]. At nuclear density they find that C a ~ 0.22 MeV -1 for one pion exchange (OPE), while it is reduced to 
C a ~ 0.08 MeV_ 1 for realistic nucleon-nucleon and N 3 LO xPT potentials. 

From Eq. (5.4) it follows that the form of the quasi-particle approximation is valid when uj <C (2it) 2 /C a (and 



uj -C Ep ) and inserting Eq. (5.4) into Eq. (5.1 ) we can obtain a low-frequency form of the zero temperature structure 
function 



A(0) 



nir l + (l + G a ) 2 (C a wy 



(5.5) 



where C a — C a /(2ir) 2 . This form of the structure function satisfies the S^ 1 sum-rule by construction but produces 
divergent results for the S® and 5*. In the following section we will combine the low and high frequency forms in 
Eq. (5.5) and Eq. (4.3), respectively, with the sum-rule constraints discussed in the preceding section to construct a 



structure function that can be used in calculations of the neutrino emissivity. 

Another commonly used limiting form of S a (u>) can be obtained by ignoring any many-particle correlations in the 
ground state. Here neutrons are distributed as free fermions in the ground state and excitations with q — and w^O 
arise as two-particle two-hole states. Following Ref. |7] we denote this as the two-body (2b) response and this is given 
by 



37m 



n 



d 3 Pi 
(2vr) 3 



(27r) 4 «5 d (P 1 + P 2 P 3 - P 4 ) S(ei + e 2 - e 3 - e 4 - w) T 2 H . 



(5.6) 



J~4 = j 1/2(1 — /3)(1 — fi), and T-L is related to the square of the matrix element defined in Eq. (5.3). For the nn 
system where only the spin-triplet two-nucleon state contributes it is explicitly given by 



n= \ \{tM' s ,v'\[S,T NN )\v,lM 3 



(5.7) 



M s Mi 



where S is the total spin and p(p') is the relative initial (final) momentum of the two-nucleon system. The matrix 
element is computed assuming plane-wave in and out states for the two-nucleon system and thus ignores many-body 
effects and initial and final state correlations. In Ref. [7] it was argued that this result can be expected to be valid up 
tou~ m w in the absence of many-body corrections. 



VI. CONSTRUCTING 

While it is clear that a unique reconstruction of S a (co) would require an infinite number of moments, here we show 
that the lowest order sum-rules and the the asymptotic forms discussed in the preceding section already provide sig- 
nificant constraints. The three sum-rules S~ x , 5°, S+ 1 provide an useful sampling of the function at low, intermediate 
and high energy respectively. Its utility in constraining the neutrino emissivity will depend on the redistribution 
of strength due to finite temperature effects. We postpone a discussion of these finite temperature effects to the 
subsequent section. Here, using as guidance the limiting forms discussed previously, we examine simple ansatze for 
the functional form of S a (u>) at T = by imposing sum-rule constraints. 

The striking feature of the sum-rule results shown in Table |T] is that at nuclear density ojq ~ Q± and is comparable 
to the Fermi energy Ep = k 2 F /2m. This suggests that the function S a {iS) contains significant strength in the region 
wo to ui%. To properly account for this we study the following simple ansatz for the frequency dependence of So-(lj) 
given by the form S'^ ow (a;) but the with a more complex behavior of t(u) given by the following forms 

t{oj) \ {uj + cj a ) 2 J \U + UJ C ) 

where the constants a,u;o and indices n,m are fit to ensure that the three sum-rules and the asymptotic forms are 
satisfied. At low frequency, this ansatz ensures that out results coincide with the results obtained by in Ref. [5] by 
Lykasov, Pethick and Schwenk where only the first term containing C a contributes. On general grounds (unitarity of 
scattering amplitudes) at large frequency, pair excitation should be quenched due to the retarded nature of nuclear 
interaction, and this quenching is naturally incorporated through the asymptotic form discussed in relation to Eq. (|4.3[) . 



8 



To better understand the sensitivity of our results to the choice of parametrization, we have also used a simple 
phenomenological form for the spin response: 

w = a a + C/^y (6 - 2) 

The high frequency tail is forced to fall off appropriately by choosing Ai — j = 9. The parameters ev,w c , and i are 
then fitted to the three sum-rules. This simple form assures that the response goes to zero at low frequency, has the 
correct high-frequency tail, and has a single peak structure. Comparisons of the two parametrizations provide some 
information on the reliability of the extracted spin response. 




Figure 2: (Color online) The spin response function S a (q = 0, u) of neutron matter at saturation density obtained by fitting 
to AFDMC sum-rules using two different ansatz are shown as the black solid and dashed curves. The inset compares the fits 
and the two-particle response at high energy obtained by confining two neutrons in a spherical cavity of radius 7 fm (red) or 8 
fm (green). The linear, low-frequency forms predicted in Ref. [36 , labeled as OPE and xPT are shown for comparison. The 
dot-dot-dashed curve is obtained using the two-body approach in Eq. (5.61 with OPE. 



Figure [2] shows the response function obtained by fitting the sum-rules and the high-energy response at saturation 
density using the two different parametrizations, Eqs. (5.5) and (6.2). For comparison, the low-frequency form of 
the structure function obtained in Ref. [53] are shown for the two choices of C, 
potentials discussed earlier. The form of the low-frequency response in Eq 



corresponding to the OPE and xPT 

In the 



(5.1) is valid only at uj -C Ep 
figure we also show the results from the two-body approach (described in Eq. (5.6)) in the Born approximation with 
OPE. At low frequency u < Ep/2, it gives similar results to the quasi-particle picture, then becomes larger at higher 
frequency since it includes the exact phase space integrals. The inset compares the fits and the two-particle response 
at high energy obtained by confining two neutrons in a spherical cavity of radius 7 fm (red) or 8 fm (green). The 
asymptotic forms and sum-rules force significantly more strength at lower energy than obtained previously. 
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The simple phenomenological fit (dashed line - Eq. (6.2 1 and the fit to the quasi-particle form (solid line - Eqs. 



(5.5) and (6.1 1) produce very similar response functions. In addition to the sum-rule constraints, we are forcing the 
response to go to zero at low frequency, have a single peak structure, and to fall off fairly rapidly at high-energy as 
obtained from the two-neutron response. Combined, these considerations place fairly tight constraints on the spin 
response of neutron matter. 




Figure 3: (Color online) The spin response function S a (q = 0, ui) of neutron matter at p = 0.12,0.16, and 0.20 fm 3 from fits 
to AFDMC sum-rules results at zero temperature . 

In Figure [3]we compare the response functions obtained over a range of densities p = 0.12,0.16 and 0.20 fm -3 . As 
expected from the sum-rules, the peak of the response functions shifts to larger energy with increasing density. The 
tensor and spin-orbit correlations are naturally of shorter range at the higher densities where the mean inter-particle 
spacing is shorter, and hence the peak shifts to higher energy. The total strength in the response is fairly flat over 
the regime of densities we consider as obtained in the sum- rule calculations for So- 

Finally, at higher density the distribution is somewhat broader as u)\ increases more rapidly with density than 
ujq. Both ujq and uj\ increase rapidly, presumably associated with the increasing importance of the shorter-range 
components of the nuclear force at and above saturation density. While we expect this trend to be qualitatively 
correct contributions due to three-body forces and from two-body currents are able to play a role in modifying this 
behavior. 



VII. EXTENSION TO FINITE TEMPERATURE AND IMPACT ON NEUTRINO PRODUCTION 

The AFDMC method we employ is restricted to zero temperature and we have not explicitly computed the tem- 
perature corrections to the sum-rules. Hence there will be several caveats to consider when using our results in finite 
temperature applications such as supernova where S a {ui) plays a role in neutrino production rates. To discus these 
we first note that there are three fundamental energy scales inherent to our present analysis of the structure function 
and the neutrino emissivity. They are: (i) typical energy at which the structure function has significant strength and 
is given by ujq and Q\\ (ii) the energy scale at which the structure function is sampled in the neutrino emissivity and 
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is denoted as lo u from the expression for Q in Eq. (2.2 1 we expect that u v ~ 5 — 6 T ; and (iii) the high energy scale 
lj c at which the asymptotic two-body behavior dominates. 

At very low temperature where w„ -C Coq and lo v < W!, the sum-rules do not provide useful constraints. Here the 
low-frequency form of S a {u)) given in Eq. (5.5) can be used to calculate the neutrino emissivity with the requirement 
that uj v <C Ep and C a uj v < 1. In practice, at nuclear density, the condition that w„ <C Ep is more restrictive and 
limits the use of the low frequency form to region where T < Ep/6. At intermediate temperature when lj v ~ loq 
or u)\ and T < Ep, the zero temperature sum-rule constraints on the form of 5° and become relevant. Here 
the temperature is intermediate and corrections to the T = sum-rules are expected to be small due to the Pauli 
principle. 

At finite temperature, the dynamic structure factor obey detailed-balance 

(7.1) 



So-(-w) = exp ( -— ) S a (u) , 



and this is reflected in Eq. (2.2) where the neutrino emissivity where the exponential term accounts for the fact 
that neutrino emission corresponds thermal fluctuations in which w is negative. There are residual temperature 
dependencies in the function S a (uj). First, from the fluctuation-dissipation theorem we have 



S a (u>) = -2 (l-cxp^))- 1 ImH». 



(7.2) 



where n R is the retarded polarization function wh ich is an odd function of to and vanishes at to = 0. To extend to 
finite temperature the zero-temperature ansatze in f VI need to be multiplied the factor (1 — exp (— lo/T))~ 1 . A second 

C a T 2 is 
We incorporate this 



source of temperature corrections arise from the fact that at low frequency the spin relaxation time T a 1 
dominated by scattering between thermally excited quasi-particles as described in Eq. (6.1) 



expected behavior by using the finite temperature expression for r CT given in Eq. (6.1 ) in the low- frequency form given 
for Sa-(uj) in Eq. (5.5). Other sources of temperature corrections exist such as those arising from transitions in which 
the excited many-particle states does not decay to the ground state and should be investigated in the future. We 
leave this for future work as it would require the development of finite temperature QMC techniques. 

With the aforementioned finite temperature extensions we employ the S ff (tjj) obtained using the sum-rule constraints 
to compute the neutrino emissivity. In Figure |4j the resulting energy loss rate Q are shown for various temperatures 
T. The large strength required by the sum-rules at intermediate energy leads to a larger neutrino emissivity compared 
to the simple extrapolation of results obtained in the quasi-particle approximation with only two-particle two-hole 
excitations. Our results are almost a factor of 2 larger than those obtained using either %PT in Ref. [36] or those 
obtained directly from nucleon-nucleon phase-shifts in Ref. [7] at T < 5 MeV. We suspect that this enhancement is 
due to correlations in the ground state that are not captured in the quasi-particle approximation. 



VIII. DISCUSSION 



Our study of the zero-temperature sum-rules suggests that the spin response function of neutron matter has a 
significant strength at energy between 40 — 60 MeV in the vicinity of nuclear density. This strength should be 
accessible at temperatures of the order of 5 — 30 MeV encountered in the supernova environment and could influence 
the rate of neutrino pair production from nucleon-nucleon processes. Although the zero-temperature sum-rules do not 
directly constrain the finite temperature response functions needed in the calculation of the neutrino emissivity they 
provide a useful guidance. For example, they can be used to test predictions obtained using quasi-particle methods 
at zero temperature, and with some caveats can be extended to low temperature where response is still dominated 
by transitions to the ground state. 

Comparisons with earlier calculations of the dynamic structure factor indicate they significantly under-predict the 
response in the regime where w ~ Ep/5 — Ep for the densities considered. There are several possible resolutions to 
this discrepancy. The quasi-particle interactions and dispersion relations used in earlier studies may not be adequate 
as finite density effects are ignored. Similarly the use of plane wave states augmented with the T-matrix, may be too 
simple to reproduce the coupling between the ground state and the excitations. Alternatively, our ansatze for the 
response function may be too simple to capture the complex structure of S a (uj). 

All of these possibilities can be studied in more detail. To quantify the interplay between increasing phase space and 
decreasing strength of the two-body interaction with increasing lo we have calculated S a {uS) in the standard approach 
using realistic potentials. Within the QMC approach there are two ways to address these issues. First, the calculation 
of higher order moments at zero temperature can provide additional constraints and test our ansatze at intermediate 
energy. Second, extensions to finite temperature will shed light on the importance of transitions not involving the 
ground state. We hope to pursue these in future work. 



11 




OPE 
sum-rule 



T [MeV] 



Figure 4: (Color online) The energy-loss rate Q at various temperature as defied in Eq. (2.2 1, for OPE, xPT and our results 
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